ApJ accepted [30 March 2004] 

Preprint typeset using I4TgX style emulateapj v. 2/19/04 



A POPULATION OF MASSIVE GLOBULAR CLUSTERS IN NGC 5128 

Paul Martini 

Harvard-Smithsonian Center for Astrophysics; 60 Garden Street, MS20; Cambridge, MA 02138; pmartini@cfa.harvard.edu 

AND 
LUIS C. Ho 

Observatories of the Carnegie Institution of Washington; 813 Santa Barbara Street; Pasadena, CA 91101; lho@ociw.edu 

ApJ accepted [30 March 2004] 



o 
o 



> 
m 
o 
o 
-t- 
o 
^1- 
o 

^ ' 
6 ■ 

a: 



ABSTRACT 

We present velocity dispersion measurements of 14 globular clusters in NGC 5128 (Centarus A) obtained 
with the MIKE echelle spectrograph on the 6.5m Magellan Clay telescope. These clusters are among the 
most luminous globular clusters in NGC 5128 and have velocity dispersions comparable to the most massive 
clusters known in the Local Group, ranging from 10-30 km s _1 . We describe in detail our cross-correlation 
measurements, as well as simulations to quantify the uncertainties. These 14 globular clusters are the brightest 
NGC 5128 globular clusters with surface photometry and structural parameters measured from the Hubble 
Space Telescope. We have used these measurements to derive masses and mass-to-light ratios for all of these 
clusters and establish that the fundamental plane relations for globular clusters extend to an order of magnitude 
higher mass than in the Local Group. The mean mass-to-light ratio for the NGC 5128 clusters is ~ 3 ± 1 , higher 
than measurements for all but the most massive Local Group clusters. These massive clusters begin to bridge 
the mass gap between the most massive star clusters and the lowest-mass galaxies. We find that the properties 
of NGC 5128 globular clusters overlap quite well with the central properties of nucleated dwarf galaxies and 
ultracompact dwarf galaxies. As six of these clusters also show evidence for extratidal light, we hypothesize 
that at least some of these massive clusters are the nuclei of tidally stripped dwarfs. 
Subject headings: galaxies: individual (NGC 5128) — galaxies: star clusters — globular clusters: general 



1. INTRODUCTION 

Globular clusters provide valuable snapshots of the forma- 
tion history of galaxies and their large sizes and luminosi- 
ties make them the most readily observable sub-galactic con- 
stituents. In addition, globular clusters exhibit surprisingly 
uniform properties that suggests a common formation mech- 
anis m. They are well-fit by single-mass, isotropic King mod- 
els (lKinglfl9r)6T) . which describe clusters in terms of scale 
radii, central surface brightness, and core velocity disper- 
sion. Detailed studies of globular clusters in our Galaxy have 
shown that in fact they only inhabit a na rrow range of the 
parameter space av ailable to King models ( Di or govskil 1 9951 
McLaughlin 2000) and other globular cluster systems in the 
Local Gro up also appear to approximately follow the same 
relations jb iorgovs ki et al.lll997t fDubath & Grillmar rjll997t 
Dubat h et alJll997t iBarmbv et alJl2002l lLarsen et al.ll2002l) . 
These parameter correlations trace a globular cluster funda- 
mental plane that is analogous to, but distin ct from, the fun- 
dame ntal plane for elliptical galaxies ( Dressier et alJl!987t 
i Diorgovski & Davislll987t iBender et al.lll992t iBurstein et all 



1997). 

Globular cluster studies that include internal kinematics 
have been confined to the Local Group due to the faint appar- 
ent magnitudes of more distant extragalactic globular clus- 
ters. These studies have therefore only included the globu- 
lar cluster systems of spiral and dwarf galaxies and not the 
globular cluster systems of large ellipticals. Yet the glob- 
ular cluster systems of ellipticals are a particularly interest- 
ing regime as they probe both a new morphological type and 
one likely to have exhibited a different and more complex 
formation history. The globular cluster systems of elliptical 
galaxies, as well as many spirals such as our own, have bi- 
modal color distributions suggestive of multiple episodes of 



formation (e.g. lKundu & Whitmorel200 It lLarsen et alJ200ll) . 
Models for the formation of these globular cluster systems 
posit that one of these populations may be the intrinsic popu- 
lation of the galaxy and subsequent mergers resulted in the 
second population either as the result of a new episode of 
globular cluster formation dS£hwei zerll987HAshman & Zepfl 
1992; Forbes. Brodie. & Grillmair 1997) or accretion of glob- 
ular cluster systems from other galaxies, including accretion 
of globular clusters by tidal stripping from other m embers of 
a cluster of galaxies (Cote. Marz ke. & Wesil998l) . 

NGC 5128 (Centarus A), as the nearest large elliptical 
galaxy, is arguably the best source for extending detailed 
globular cluster studies outside of the Local Group. While 
NGC 5128 is the central galaxy of a large group, rather than 
a giant elliptical in a cluster, it likely had a similar formation 
history to its larger cousins. The most relevant similarity is 
the strong evidence for a recent, gas-r ich merger (for a recent 
review of NGC 5128 see Israeli Il998l) . Estimates of the size 
of the NGC 5128 globular cluster population suggest that it 
has a total of ~ 2000 clusters, approximat ely a fact or of 3 
more than the entire Local Group (Harris et al. 1984). Sim- 
ple scaling arguments suggest that NGC 5128 should possess 
a number of extremely massive globular clusters and there- 
fore is not only a good target for study of the globular clus- 
ter system of an elliptical galaxy, but also for study of how 
well the fundamental plane relations established locally apply 
to more massive globular clusters. A recent p hotomet ric and 
spectrosc opic study of NGC 5128 by Peng. Ford. & Fr eeman! 
(2004a b) concluded that the metal-rich globular clusters may 
have a mean age of Gyr, while an analysis of their pho- 
tometric data y ields a metallicity range of -2.0 through +0.3 
dYi et al J2004I) 

The most massive globular clusters can also be used to ex- 
plore connections between the formation processes for star 
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clusters and galaxie s. While fundamental plane studies (e.g. 
iBurstein et al.i n"997) clearly illustrate a significant mass gap 
between the most massive Galactic globular clusters and the 
least massive dwarf galaxies, there have been encroachments 
into this gap from both sides. For many years, studies have 
speculated that at least some globular cluster s may be the re- 
mains of tidally- stripped dwarf galaxy nucl ei JZinnecker et alJ 
119881 lFreemanlfT993t iBassino et alii 1994) . Two of the most 
massive globular clusters in the Local Group, wCen in our 
Galaxy and Gl in M31, h ave been interpreted as the nuclei 
of tidally stripped dwarfs jMevlan et alJl200"Tl iGnedin et alJ 
120021 iBekki & Freemanl 120031) . From the galaxy side, re- 
cent studies of nucleated dwarf galaxies in th e Virgo cluster 
JGeha. Guhathakurta. & van der Marel l2002j) and ultracom- 
pact d warf galaxies in the Fornax cluster fDrinkwat er et alJ 
2003) show some similarities between these least-massive 
galaxies and the most massive globular clusters. This mass 
gap may thus reflect the scarcity of the most massive globu- 
lar clusters and the difficulty of kinematic measurements for 
the least massive, lowest surface-brightness dwarf galaxies, 
rather than a physical separation. 

In this paper we present velocity dispersion measurements 
for 14 globular clusters in NGC 5128. These globular clus- 
ters we re selected from th e Hubble Space Telescope (HST) 
study of H arris et al]d2002) and therefore have well-measured 
structural parameters. These data are combined to estimate 
masses for these clusters, masses that are among the largest 
known for any star clusters and comparable to the nuclei of 
the lowest-mass galaxies. In the next sections we describe the 
observations, data processing, and velocity dispersion mea- 
surements. Analysis of these measurements is described in §5 
and the potential link between star clusters and galaxies is ex- 
plored in §6. Our results are summarized in the final section. 
Throughout this paper we adopt th e distance of 3. 84 ± 0.35 
Mpc for NGC 5128 determined by iReikubal (120041) from the 
brightness of the tip of the red giant branch and the Mira 
period-luminosity relation. 

2. OBSERVATIONS 

Spectra of 14 globular clusters in NGC 5128 were ob- 
tained in the course of seven nights in March 2003 with the 
MIKE echelle spectrograph ( Bernstein et alJ 120031) and the 
6.5m Magellan Clay telescope at Las Campanas Observatory 
(see Table [0. MIKE is a double echelle spectrograph capa- 
ble of resolution R = 19,000 on the red side when used with 
the 1" slit selected for these observations. This corresponds 
to a velocity resolution FWHM of 15.8 km s" 1 or a = 6.2 km 
s -1 . MIKE was located at the East Nasmyth port but not di- 
rectly attached to the instrument rotator. The advantage of 
this configuration is that the instrument is gravity invariant 
and therefore potential calibration difficulties due to instru- 
ment flexure are completely avoided. The disadvantage of 
this configuration is that targets could not be observed at the 
parallactic angle and that they rotate with respect to the spec- 
trograph slit. In any event, this configuration was the only 
one available at the time. To minimize lost light due to at- 
mospheric dispersion, the position angle of the spectrograph 
was set to correspond to the parallactic angle for an object at 
airmass 1.3; below airmass ~ 1.4 the atmospheric dispersion 
correction is effectively negligible and most of our observa- 
tions took place in this airmass range. Any lost light due to 
atmospheric dispersion is unlikely to significantly impact our 
velocity dispersion measurements as the velocity dispersion 
gradient of the clusters is not large. The rotation of the slit 



TABLE 1 

Globular Cluster Observations and Measurements 



ID 


V [mag] 


Exptime [s] 


SNR 


n -1 i 

<t [km s J 


C2 


18 33 


10800 


[ j 


14 2 ± 2.0 


C7 


17 10 


5400 


16 


22 4 ± 2 1 


Cll 


17.70 


9000 


18 


17.7 ±1.9 


C17 


17.61 


7200 


19 


18.9±2.0 


C21 


17.77 


10800 


12 


20.8 ±1.9 


C22 


18.14 


10800 


15 


19.1 ±2.0 


C23 


17.19 


7200 


25 


31.4±2.6 


C25 


18.33 


12600 


5 


12.2 ±2.0 


C29 


17.94 


14400 


12 


16.1 ±2.1 


C31 


18.37 


12600 


8 


15.0±2.1 


C32 


18.31 


10800 


7 


15.8 ± 1.4 


C37 


18.34 


10800 


10 


13.5 ±1.6 


C41 


18.56 


19800 


8 


9.6 ±2.0 


C44 


18.61 


18000 


8 


9.1 ±2.0 



NOTE. — Present sample of NGC 5 128 globular clust ers. Cols. 1 
and 2 li st the cluster ID and apparent V magnitude from Pe ng et all 
1 2004a) and col. 3 our exposure times. The mean signal-to-noise 
ratio over the range 5000 - 6800 A per 0.2 A is listed in col. 4. The 
measured velocity dispersions and uncertainties, described in ^4] are 
listed in column 5. 



also has a negligibl e impact as all of these globular clusters 
are fairly round (Harr is et all2 002). 



While MIKE was used to obtain both red and blue data, suf- 
ficient signal-to-noise ratio (SNR) observations of these (in- 
trinsically red) clusters were only obtained on the red side. To 
enhance the SNR, the observations were binned 2x2, leading 
to 0.286" pixel" 1 and a velocity scale of 4.2 km s _1 pixel _I . 
Table H] lists al l 14 c lusters according to their identification 
in lHarris et al J (120 021). a long with their apparent V magnitude 
from lPeng et alJ ( 2004a), the total integration time, and the av- 
erage SNR per pixel over the wavelength range 5000 - 6800 
A when binned to 0.2 A pixel" 1 . The total integration time 
listed is the sum of an integer number of 1800 s exposures. 
This exposure time was chosen to facilitate cosmic ray re- 
moval, yet also avoid a significant contribution from detector 
noise. Spectra of 16 velocity dispersion templates were also 
obtained, ranging in spectral type from G through M stars and 
luminosity classes from IV through I. These templates were 
observed in a similar manner to the globular clusters. Typi- 
cally three integrations of a few seconds each were obtained 
for each template. Arc lamps for wavelength calibration were 
obtained before or after the observations of each target. Sev- 
eral sets of flatfield exposures using a quartz lamp were also 
obtained over the course of the seven-night observing run. 

3. DATA PROCESSING 

The data were processed with a combination of tasks from 
the IRAF 1 echelle pac kage a nd the sky-subtraction software 
described by Kelson) ll2003l) . All images were overscan- 
subtracted and a simple cosmic -ray cleaning algorithm was 
applied to the individual cluster exposures. The images were 

1 IRAF is distributed by the National Optical Astronomy Observatories, 
which are operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Science Foun- 
dation. 
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also rotated such that the dispersion direction was approxi- 
mately parallel to the rows of the detector. Inspection of the 
several sets of flatfield exposures obtained over the course of 
the run showed no measurable variation and all of the flatfield 
exposures were therefore summed to produce one master flat- 
field. This flat was then normalized with low-order fits to each 
order and divided into all of the object exposures prior to sky 
subtraction. 

Sky subtraction of MIKE data is complicated by the nonlin- 
ear transformation between the orthogonal (dispersion, spa- 
tial) and "natural" detector (row, column) coordinate systems 
and the slight undersampling of the data. The marginally un- 
dersampled sky lines are tilted with respect to the detector 
coordinate system and this tilt angle vari es as a functio n of 
spatial position. As discussed in detail by Kelson (2003|), the 
classical approach of rectifying and rebinning the image prior 
to sky subtraction can result in significant residuals in the case 
where the sky lines are critically- or undersampled. Superior 
sky subtraction can be obtained if a sky model is calculated 
and subtracted prior to any rectification and rebinning of the 
data. 

The key to this technique is to calculate the transformation 
between the array coordinates (x,y) and the orthogonal disper- 
sion and spatial coordinate system (x r ,y t ). By collapsing the 
spectrum along the spatial coordinate y t , the profile of each 
sky line can then be measured as a function of the coordi- 
nate x r in which the sky line is well-sampled (due to the tilt 
of the observed spectrum). A two-dimensional sky model can 
then be constructed with the same pixelization as the origi- 
nal data, thus avoiding the residuals that result from rebin- 
ning marginally or undersampled data. To employ Kelson's 
method, we first calculated the transformation y, = Y(x,y) with 
his getrect command, using a flatfield frame to trace the 
boundaries of each order. These boundaries were defined with 
f indslits (treating them as individual slits in a multislit 
mask). We then calculated x r =X(x,y t ) with getrect ap- 
plied to each order. Accurate calculation of this transforma- 
tion requires several lines per order. As there are few airglow 
lines in the higher orders (shorter wavelengths), we created a 
composite image of spectroscopic lines by combining a very 
high SNR sky frame (calculated by summing a large number 
of our exposures of the fainter clusters) with ThAr and NeAr 
lamp spectra. As many ThAr lines saturated in the 13 reddest 
orders, these orders were masked out before construction of 
the composite line image for getrect. A sky model image 
for each globular cluster exposure was then calculated with 
skyrect and subtracted from the flatfielded data. Observa- 
tions of the template stars were sufficiently short that no sky 
subtraction was deemed necessary. 

The dispersion solution was calculated with the 
ecidentify task in the IRAF echelle package and 
ThAr lamp spectra. Although most of the reddest orders were 
saturated by bright lines as described above, a combination 
of multiple lines in the bluer orders and several unsaturated 
lines in the redder orders led to a good solution (rms ~ 0.005 
- 0.007 A) with a fourth-order Legendre polynomial in x 
and y. All of the object spectra were then extracted with 
apall. Because most of the individual cluster spectra were 
too faint for a good trace, particularly in the bluest orders, 
we used a reference trace from a brighter object. This trace 
was recentered based on the centroid of the cluster trace in 
the reddest, highest SNR orders. Each individual exposure 
was extracted in this manner as several-pixel shifts were 



observed between a small number of the multiple exposures 
on individual clusters. The width of the spectral extraction 
window was 9 pixels or 2.6". The arc lamp was extracted 
for each object with the same trace and the dispersion 
solution was calculated with the ecreidentify task. The 
spectra were then placed on a linear wavelength scale with 
the dispcor task for measurement of their line-of-sight 
velocity dispersion a. Template stars and the flux standard 
HR 1544 were extracted in a similar manner. Orders 68 - 
39 of C23, the cluster with the highest SNR, are shown in 
Figure ^ The calcium triplet region (order 40) is shown 
for all 14 clusters in Figure |2] No correction for telluric 
absorption was applied to these data. 

4. VELOCITY DISPERSION MEASUREMENT 

The two most commonly employed techniques for veloc- 
ity dispersion measurement are direct-fitting in pixel space 
(e.g. lRix & WhiteH992tlKelson et al.l2000tlBarth et alJ2002t) 
and cross-correl ation techniques in the Fourier domain (e.g. 
Illing worth! 1 1 9761) . After experiments with both techniques 
we chose the cross-c orrelation method developed in detail by 
Tonrv & Davisl (11 9791) as implemented by the rvsao pack- 
age (Kurtz & Mink 1998) in IRAF. The main advantage of 
the cross-correlation function (hereafter CCF) is that it is less 
sensitive to the relative line strengths in the targets and tem- 
plates. This is likely to be the case for the present study due 
to metallicity differences between the globular clusters and 
the Galactic template stars observed. Trial implementations 
of direct-fitting in pixel space confirmed that most of the stel- 
lar templates are not good matches in detail to these globular 
clusters. Direct pixel-fitting did show that the observed stel- 
lar templates with spectral types between G8 and K3 were the 
best matches to the globular cluster spectral energy distribu- 
tions. The stellar templates used for this study are: HD 80499 
(G8III), HD 95272 (K0III), HD 88284 (K0III), HD 43827 
(K1III), HD 92588 (K1IV), HD 44951 (K3III), and HD 46184 
(K3III). 

The first step in implementation of the CCF technique is the 
calculation of the relation between the CCF width and the ob- 
served line-of-sight velocity dispersion a. To measure this re- 
lationship, we convolved a range of the observed stellar tem- 
plates by Gaussian functions with a = [5,7,10,14,20,28,40,56] 
km s"'and used a simple spline fit to calculate the relation be- 
tween CCF width and input a. We calculated this relationship 
for each of the spectroscopic orders and determined that the 
uncertainty in a measurements in orders with significant tel- 
luric absorption (50, 48, 47, 45, 42, and 41) was considerably 
higher than those without. Given the large number of orders 
available, these orders were excluded from subsequent anal- 
ysis. Based on measurements from the remaining orders we 
determined that all of the globular clusters have a in the ap- 
proximate range of 10 - 30 km s" 1 . The velocity dispersions 
of these clusters are thus well-resolved by this instrument con- 
figuration. 

While the CCF technique is known to be relatively insen- 
sitive to template mismatch, we tested the importance of this 
potential source of systematic error by convolving each spec- 
tral type by the range of velocity dispersions given above and 
using other templates of similar but not identical temperature 
to measure a (e.g. a K0III template to measure an artificially 
broadened K3III template). This analysis showed that tem- 
plate mismatch produced at most an rms uncertainty of 1 km 
s" 1 when a was calculated for a single order. Given the large 
number of orders available and our estimate of the range in a 
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Fig. 1. — Orders 39 - 68 of the NGC 5128 globular cluster C23. Each order has been divided by a low-order polynomial fit to normalize 
the spectrum (dark solid line, left axis). The SNR as a function of wavelength is shown below the spectrum (light solid line, right axis). The 
smooth SNR variation is due to the transmission of the individual orders, which have been stiched together. 
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Fig. 2. — Calcium triplet region (order 40) for all 14 NGC 5128 globular clusters. 



of these clusters, we then experimented with combining mul- 
tiple orders together as well as binning the data to improve 
SNR. A range of tests showed we obtain excellent and very 
stable results over the wavelength range 5000 - 6800 A and 
the spectra rebinned to a scale of 0.2 A pixel" 1 . The rms vari- 
ation due to template mismatch over this range was less than 
0.5 km s . CCFs for all 14 clusters are shown in Figure [5] 
and the measured a and estimated uncertainties are listed in 
Tabled 

The 5000 - 6800 A spectral range does not include the 



higher SNR CaT region. While the CaT lines are quite strong 
for all of the clusters (see Figure |2J, these lines are in fact 
so strong that they exhibit significant damping wings. These 
damping wings broaden the line profiles of the template stars 
significantly and they are not a good representation of the in- 
strumental line profile. The CaT lines in the template stars 
have intrinsic widths comparable to or larger than the mea- 
sured velocity dispersions of the globular clusters over the 
5000 - 6800 A region. We therefore did not use these lines 
to measure the cluster velocity dispersions. The neighboring, 
red orders have comparable SNR to the CaT region, although 
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these orders are contaminated by significant telluric absorp- 
tion. Although the telluric absorption is to some extent cor- 
rectable, these orders were also discarded for the velocity dis- 
persion measurement because the SNR in the 5000 — 6800 A 
region is sufficiently high that SNR is at most a modest con- 
tributor to the total error budget (see below). While there are 
still reasonably strong lines in the 5000 - 6800 A region, no- 
tably Mgb and Ha, these lines only make a minor contribution 
to the velocity dispersion measurement. 

Given the range of SNR and a of these data, we performed 
simulations with our template stars to determine the quality 
of our measurements as a function of both SNR and a. These 
tests were performed by first convolving the best-fitting tem- 
plate (HD 88284) with the range of velocity dispersions de- 
scribed above. We then added noise to these convolved spec- 
tra to produce a set of output spectra with average SNR of 
[2,4,8,16,32,64,128] for each velocity dispersion. As the SNR 
per pixel is quite variable due to both the blaze function and 



the red color of these clusters, we used an input SNR spectrum 
from one of the cluster observations (C23, shown in Figure[l) 
to reproduce the pixel-to-pixel variation in the SNR, while 
forcing the mean SNR in the spectra to equal the specified 
values. The widths of these input spectra were then calculated 
in a similar fashion to the globular cluster measurements de- 
scribed above. The results of this analysis for a = 10-28 km 
s (the range most appropriate to this study) are shown in Fig- 
ure[4] This figure displays the difference between the input a 
and the measured a at a given SNR, normalized by the input 
a. 

As expected, the deviation between the input a and the 
measured a increases toward lower SNR. In particular, these 
simulations show that at lower SNR the measurements will 
tend to overestimate a. This may be because the addition 
of noise systematically erases narrower features more read- 
ily than broader features. The simulations also show that at a 
given SNR the velocity dispersion measurement will be more 
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Fig. 4. — Velocity dispersion residuals as a function of signal-to- 
noise ratio for the spectral region 5000 - 6800 A and binned to 0.2 
A pix~'. The a values shown closely correspond to the observed 
range for the NGC 5128 globular clusters. Less than half of the 
cluster observations have mean SNRs below ten. Tickmarks along 
the dashed line mark the location of the model SNR values. Only 
the panel for a = 10 km s plots a residual value for SNR= 2. The 
residuals were larger than 30% at SNR= 2 for the remaining three 
velocity dispersions. 



accurate if the intrinsic velocity dispersion is lower. This is 
likely a reflection of the net SNR over the number of pixels 
that comprise a given feature. The fewer pixels that span, e.g. 
the absorption lines employed here, the higher the net SNR of 
the feature relative to a spectrum of the same mean SNR per 
pixel with a larger a. This systematic effect on the a measure- 
ment impacts our observations at most at the 5% level. Three 
of the clusters with velocity dispersions of on order 10 km 
s~'(C25, C41, C44) have a mean SNR less than 10 over then 
5000 - 6800 A range in our rebinned spectra, while three (C2, 
C31, C32) in the 14 km s"'range have a mean SNR ~ 10 (see 
Table [3. The measured velocity dispersions for these low- 
SNR clusters have been adjusted to account for the expected 
overestimate. The remaining clusters do not suffer from sys- 
tematic errors due to low SNR. 

The resolution of our rebinned spectra could also impact the 
measured a as the velocity dispersion resolution now varies 
from 12 km s"'at 5000 A to 8.8 km s _1 at 6800 A. How- 
ever, measurements of a with a dispersion of 0.1 A pixel" 1 
yield consistent results with the values in Table [0 although 
with somewhat larger scatter due to the lower SNR. The re- 
lationship between measured a and SNR at each a shown in 
Figure 0] also demonstrates that the measured a converges to 
the input value as the SNR increases. The measured values 
quoted in Table[^also agree well with measurements obtained 
from single, unbinned orders and direct-pixel fitting, although 
the latter approaches have larger scatter. We note that prelimi- 
nary velocity dispersions measurements from high-resolution 
spectroscopy of four NGC 5128 globu lar clusters w ere ob- 
tained by Dubath and reported in Harri s et all (|2()02). These 
measurements are marginally consistent with our own. 

The 1" x 2.6" spectral extraction window we employed cor- 
responds to a physical area of 18.6 x 48.4 pc, which is larger 



than the half-light radius for all of these clusters (e.g. see Ta- 
ble|2j. In King models for globular clusters, the velocity dis- 
persion is a function of radius. Our measurements, because 
they include most of the cluster, should correspond well to the 
global, one-dimensional velocity dispersion a^. In contrast, 
the core one-dimensional velocity dispersion ctq is approxi- 
mately the quantity measured for Galactic globular clusters 
and the most commonly quoted quantity in the literature. The 
core value 00 is larger than in the King models and the 
degree of difference depends on cluster concentration. For 
the clusters studied here, this correction is on order 5 - 10% 
for the range of cluster concentrations in this sample (see e.g. 
Diorgovs ki et aljfl997l) . The errorbars quoted in Table [T] in- 
clude the uncertainties in this correction, which are caused by 
un certainties in the m easured structural parameters reported 
bv lHarris et alJ d2002l) and the exact placement of the slit, as 
well as the uncertainties due to template mismatch and SNR. 

5. ANALYSIS 

5.1. Observed Cluster Properties 

These 14 globular clusters are the brightest clusters with 
King model struct ural parameters measured with HST by 
Harris et alJ ll2002l) . These models are described by the pa- 
rameters Wo (central potential), r c (core radius) and c (clus- 
ter con centration, where c = \ogr,/r c and r, is the tidal ra- 
dius). lHarris et"al1 ll2002l) also fit for ellipticity. Both r c 
and r/, (the model half-mass radius) from their fits are pro- 
vided in Table |2j converted to parsecs, and their measure- 
ment of c is also given. The projected half-light radius is 
Rh w 0.73 r/,. As these values span the approximate range 
of 1.4 < c < 2.1, none are core-collapse clusters. All 14 
of these clusters w ere also included in the recent study by 
IPeng et alJ J2004al) . who have obtained UBVRI photometry 
and radial velocities for over 200 NGC 5128 globular clus- 
ters. As most of the Harris et al. (2002) photome try is from 
unfiltered STIS images, we have listed instead the lPeng et alJ 
(2004a) Vo, (B-V)o, and (V-I)q values in the table, where 
these quantities have been dereddened based on their quoted 
E(B — V) and assuming Ry = 3. 1. We also include the cen- 
tral surface brightness fiy from Harris et al. (2002), but re- 
calibrated by the differe nce between their quoted V magni- 
tudes and those in P eng et all ll2004al) . and have calculated the 
mean, reddening-corrected surface brightness within the pro- 
jected half-light radius: {nv)h = Vb + 0.75 -2.5 log -kR\. All 
of these dereddened values only account for Galactic dust and 
do not include any dust intrinsic to NGC 5128. While most of 
these clusters are not near NGC 5128's prominent dust lane, 
dust may affect the colors of some of these clusters. 



The co lor distribution of NGC 5128 glob ular clusters is 
bimodal (iHeld et alJl 19971 IPeng et aljl2004al) . s imilar to the 
globular cluster systems of many galaxies (e.g. lLarsen etall 
1200 11 iKundu & Whitmorell2001l) . A bimodal color distribu- 
tion is most commonly ascribed to a bimodal distribution in 
metallicity as color is a reasonable proxy for metallicity in 
old single stellar populations. In the case of NGC 5128, this 
galaxy appears to have undergone a major merge r in the recent 
past. Based on the strength of H/3, IPeng et all l l2004al) con- 
clude that the more metal-rich glo bular clusters formed 5^ 
Gyr ago, while[Kavirai et al. (2004) derive an age estimate of 
1-2 Gyr for the metal-rich population in their s tudy of the 
metallicity distribution function (derived from the lPeng et alJ 
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TABLE 2 

Structural and Photometric Measurements 



ID r c [pc] r„[pc] c $ {w)h V M v (B-V) (V-I)o 



C2 0.80 8.49 1.99 16.27 17.70 18.10 -9.82 0.70 0.86 

CT 1.41 10.00 1.83 15.46 16.79 16.83 -11.09 0.75 0.91 

Cll 1.30 10.41 1.88 16.19 17.58 17.54 -10.38 0.93 1.11 

C17 2.27 7.60 1.43 16.14 16.65 17.29 -10.63 0.77 0.88 

C21 1.21 9.27 1.86 15.97 17.32 17.53 -10.39 0.78 0.93 

C22 1.10 5.06 1.62 15.60 16.29 17.81 -10.11 0.79 0.91 

C23* 0.87 4.41 1.67 14.37 15.06 16.88 -11.04 0.96 1.10 

C25* 0.99 8.32 1.90 16.34 17.71 18.15 -9.77 0.96 1.14 

C29 1 ' 1.19 9.16 1.87 16.21 17.52 17.75 -10.17 0.87 1.05 

C31 0.91 5.29 1.74 15.71 16.58 18.01 -9.91 0.92 1.12 

C32* 0.56 7.24 2.06 15.88 17.33 18.07 -9.85 0.94 1.10 

C37 1 0.58 4.43 1.87 15.31 16.28 18.09 -9.83 0.84 0.99 

C41 0.78 5.99 1.87 15.94 17.04 18.19 -9.73 0.87 1.06 

C44 1.25 7.60 1.70 16.45 17.68 18.32 -9.60 0.68 0.84 



NOTE. — Structural and photometric parameters of our sample based on measurements from the 
literature. For each cluster in column 1 we list r c , ri,, c, and [i v from Harris et al. 1 2002) in cols. 2-5, where 
we have adopted a distance to NGC 5 128 of 3.84 Mpc. Col. 6 contains the mean surface brightness within 
the projected half-light radius, while cols 7-10 are the reddening-corrected Vo magnitude, the absolute V 
magnitude, and the (B — V)q and (V — I)n colors from Pens et al. 1 2004a). We have used this photometry 
to apply a c orrection to /x?, (see text) and calculate (jj,v)h- The six clusters marked with an x superscript 
indicate that Harris et al. 1 2002) found evidence for extratidal light in their surface photometry. 




0.4 0.6 0.8 
(B-V)„ 

Fig. 5. — Reddening-correcte d colors of NGC 5 128 globular clus- 
ters. All clusters observed by Peng et al. 1 2004a ) are shown in the 
histogram, while the globular clusters in the present study are in the 
hatched histogram. The latter clusters span a range in color repre- 
sentative of the NGC 5128 globular cluster system. 



J2004al) U and B photometry). Figure |5] shows the (B-V)o 
and (V - I)o color distrib utions f or all of the globular clus- 
ters studied by Peng et al. (2004a). The clusters in the present 
study are represented by the hatched histogram and span the 
full range in color of the majority of the NGC 5128 globular 
cluster system. Figure |6] demonstrates that these clusters do 
not show any strong correlations between color and velocity 
dispersion, central concentration, mass, or mass-to-light ratio 
(see next section). 

5.2. Cluster Masses and Mass-to-Light Ratios 

One of the most valuable and readily calculated properties 
of a globular cluster is its total mass. There are multiple ways 
to estimate the mass of a globular cluster, although the most 
straightforward approach is to use the virial theorem: 



:2.5- 



• n, 



(1) 



where (v 2 ) ss 3tr 2 is the mean-squa re speed of the c luster stars 
and r/, is the half-mass radius (Bin nev & Tremainell987l) . The 



TABLE 3 

Masses and Mass-to-Light Ratios 



ID 


M vir [Mq, 6 ] 


M K mg [M 0j6 ] 


fJ-King 


TV© log 


Eb [erg] 


C2 


2.99 ±0.67 


1.74 


52.4 


4.1 ±1.6 


52.16 


CT 


8.75 ±1.46 


5.85 


39.7 


3.7 ±0.9 


53.10 


Cll 


5.69 ±1.04 


3.64 


43.0 


4.7 ±1.3 


52.68 


C17 


4.73 ±0.86 


3.84 


21.7 


3.1 ±0.7 


52.82 


C21 


7.00 ±1.15 


4.52 


41.7 


5.7 ±1.5 


52.92 


C22 


3.22 ±0.59 


2.43 


28.6 


3.4±0.8 


52.61 


C23 v 


7.59±1.21 


5.61 


30.9 


3.4 ±0.7 


53.40 


C25- r 


2.16±0.52 


1.35 


44.5 


3.1 ± 1.2 


51.93 


C29 r 


4. 14 ±0.87 


2.71 


42.3 


4.1 ±1.3 


52.47 


C31 


2.07 ± 0.47 


1.48 


34.4 


2.6±0.8 


52.16 


C32- r 


3.15±0.51 


1.74 


60.8 


4.2 ± 1.2 


52.24 


C37* 


1.41 ±0.28 


0.92 


42.3 


1.9±0.6 


51.85 


C41 


0.96 ±0.30 


0.63 


42.3 


1.4 ±0.7 


51.39 


C44 


1.10±0.36 


0.70 


32.3 


1.8±0.9 


51.41 



NOTE. — Derived parameters for NGC 5128 globular clusters. For 
each cluster in column 1, the virial mass and King model masses (in units 
of 10 6 Mq) are listed in columns 2 and 3. The dimensionless King mass, 
mass-to-light ratio, and binding energy are provided in columns 4-6. 



virial masses for these clusters (in units of 10 6 M©) are listed 
in Tabled 



An alternative approach is to use the King para meters to 
calculate the cluster mass. From Illingw ortr] i 1 9761) . the mass 
is 



Mi 



King 



= 167 r c yttcr , 



(2) 



where fi is 
1966) and 



the 

£Q_' 



dimensionless King model mass (King 
the core velocity dispersion (see also 
Richst one & Tremainelll9 86). As discussed previously, we 



have e stimated an from the cluster concentration parameter 
(Diorg ovski et aljfl997l) . The dimensionless mass /i is also a 
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Fig. 6. — Reddening-corrected V-I color of NGC 5128 globular clusters as a function of velocity dispersion (upper left), central concentration 
(lower left), virial mass (upper right), and mass-to-light ratio (lower right). All NGC 5128 clusters are marked with circles, while clusters 
identified bv lHarris et alj {2002) to possibly have extra tidal light are marke d with partly filled circles. The dashed line in the lower-left panel 
represents the mean core mass-to-light ratio derived by McLaughlin (2000) for Galactic globular clusters. The 0.18 magnitude uncertainty in 
the NGC 5128 distance modulus iReikuba 2004) translates into an additional, systematic uncertainty of ~ 18% in the mass-to-light ratio. wCen 
(cross) and Gl (star) are also shown in this panel. 



function of concentration. The King model masses for these 
clusters are ~ 50% less than the virial mass estimates. It 
is not unusual for these two estimates to disagree (e.g. see 
iMevlan et all200il for wCen and Gl). The uncertainties (not 
shown) in the King masses are larger than the virial mass es- 
timates due to uncertainties in the measurement of c and r c , in 
addition to the uncertainties in a. We therefore use the virial 
mass in the subsequent discussion. 

Once the mass has been estimated, the V-band mass-to- 
light ratio in solar units Yy© is readily calculated. The virial 
mass is divided by the V-band luminosity Ly iQ , which is com- 
puted from the reddening-corrected magnitude Vo (assuming 
Mv,q = 4.83). The mass-to-light ratios for these clusters are 
listed in column 5 of Table[3] These values fall in the range 1 .4 
- 5.7 and have an average value of (Yy ©) ~ 3. This average 



is larger than the mean core mass-to-light ra tio of 1 .45 ± 0.1 
for globular clusters in the Galaxy ( McLaughlin 1200(J) and 
the comparable core and glob al values of 1.53 ± 0.18 from 
four globular clusters in M33 (Larsen et al. 2002). If we had 
adopted the King mass estimates to compute the mass-to-light 
ratios, rather than the virial mass estimates, the mass-to-light 
ratios would be ~ 50% less. The difference between the 
NGC 5 128 globular clusters and those in the Local Group may 
be due to either intrinsic differences or due to measurement 
uncertainties. If the latter is the case, then the masses are over- 
estimated, the luminosities are underestimated, or both. An 
overestimate in the mass could be explained by either an over- 
estimate of a by ~ 50% or an overestimate of r/, by a factor of 
2, although both of these possibilities vastly exceed the esti- 
mated uncertainties in these quantities. The mass-to-light ra- 
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tios between the different globular cluster systems could also 
be brought into agreement if the luminosities are underesti- 
mated by a factor of 2. This exceeds the estim ated uncertainty 
in the distance to NGC 5128 jReikubal20 04). although intrin- 
sic reddening could also contribute and NGC 5128 is well- 
known for its prominent, large dust lanes. The luminosities 
could also be underestimated if there is a significant contribu- 
tion to the surface brightness at large radii, where the SNR is 
poor. While a conspiracy of all of these potential sources of 
error could bring the mean mass-to-light ratios of these glob- 
ular cluster systems into agreement, we conclude that the dif- 
ferences are likely real. 

Further support of these large mass-to-light ratios is pro- 
vided by some of the most massive globular clusters in the 
Local Group, including the Galactic globular cluster uiCen, 
which has T v ^ = 2 4 - 3.5 and Gl in M31 with T = 3.6-7.5 
(Mevl an et alJ 119951 1200 ll) . The lower-right panel of Fig- 
ure [6] shows the mass-to-light ratios for the NGC 5128 glob- 
ular clusters as a function of (V - 7)q. No obvious trend with 
color is apparent, as might be expecte d if the generally re dder, 
more m etal-rich clusters identified by Peng, Ford, & Freeman 
(2004b) formed more recently. The massive, Local Group 
clusters usCen and Gl are also shown for comparison. For 
these cluster s the average of the virial and King mass es- 
timates from M evlan et alJ ( 120011) were used for the mass- 
to-l ight ra t io, w hile the integrated colors ar e either from 
the lHarrisI ll 19961) compilation (wCen) or from Heaslev et al.l 
(119881) . 

5.3. Fundamental Plane 

NGC 5 128 has a significantly larger population of the most 
massive and luminous globular clusters than any galaxy in the 
Local Group. Our cluster sample is dominated by these clus- 
ters due to observational constraints, yet this population is ar- 
guably the most interesting as they provide an opportunity to 
verify and extend relationships for local clusters to more ex- 
treme examples of the population. The clusters in our sample 
are comparable in mass to wCen, th e most massive G alactic 
globular cluster at M = 5 x 10 6 M^ dMevlan et alJll995l) and 
M31' s Gl, which with M = (7 - 17) x 10" M^ JMevlanetaTI 
120011) is the most massive cluster known in the Local Group. 

The main relationship of interest is the globular cluster fun- 
damental plane, the approximately two-dimensional structure 
occupied by clusters in the three-dimensional space defined 
by, e.g. cent ral velocity dispersion, surface brightness, and 
core radius. Diorgovski ( 1995) demonstrated that the plane 
occupied by globular clusters is consistent with the expecta- 
tions for virialized cluster cores, 

r c ~a 2 XT-\ (3) 

and a constant mass-to-light ratio. This appears to be the case 
for both the core properties as well as properties derived at 
the half-light radius, with the differences between these two 
regimes due to the degree of central concentration. 

In the top panels of Figure [7] w e plot two projections of 
the core fundamental plane from Diorgovski (1995) for the 
NGC 5128 globular clusters, along with clusters in the Milky 
Way, M31, M33, and the Magellanic Clouds. These projec- 
tions show that the NGC 5128 clusters are relatively large 
and have higher surface brightnesses than Galactic clusters, 
although are more similar to globular clusters studied in other 
Local Group galaxies. Half-light projections {bottom pan- 
els) show all of these globular cluster systems follow similar 



trends to the Milky Way system, although extragalactic globu- 
lar clusters appear to have systematically lower mean surface 
brightness within the half-light radius for a given a. This may 
reflect a bias toward selection of bright objects that neverthe- 
less appear marginally resolved in ground-based images or a 
relative overestimate of R/, for more distant clusters. 

Numerous literature sources were employed to obtain data 
for globular clusters in the L ocal Group. Data for the 
Milky Way were obtained from Pryor & Mevlan ( 1993) and 
Harris (19 9g). Velocity dispersion da ta f or M31 were ob- 
tained from lDubath & Grillmairl i ll 9971) and lDiorgovski et all 
119971), while struc tural parameters were obtained from 
Barmbv"et al.l (120 021) or earlier wo rk ( B atti sjirri£tayjl982l 
Crampton et al. 1985; Bendinelli et al. 1993; Fusi Pecci et al. 
Il994t iRich et alJll996l iGrillmair et alJl!996l). Data for the 
Magellanic Clouds were obtained from |Djib_ath_g£ al J (119971) 
and data for M33 were obtained from Larsen et al. (2002), 
while the photometry for the latter was dereddened with the 
E(V - I) values derived by ISaraiedini et al.l (119981) . Data 
for the three Fornax dwarf galaxy's globular clusters with 
velocity dispersion measurements and structural para meters 
were obtained from iDubath. Mevlan. & Mavorl (119921) and 
Mac kev & Gilmorel (120031) . 

McLaughlin (2000) has recently recast discussion of the 
fundamental plane for globular clusters in terms of mass- 
to-light ratio, luminosity, binding energy, and Galactocen- 
tric radius. In the context of the four independent param- 
eters that characterize single-mass, isotropic King models, 
he found that Galactic globular clusters could be described 
with a constant mass-to-light ratio and a binding energy reg- 
ulated only by luminosity and Galactocentric radius: E/, = 
7.2 x 10 39 (L/L Q ) 2 05 (r yc /8kpc)- ' 4 erg. Figure |8] shows that 
NGC 5128 clusters follow approximately the same relation 
between binding energy and luminosity as Milky Way and 
other globular clusters in the Local Group, although the 
NGC 5128 and other clusters with higher estimated M/L fall 
above the relation traced by Milky Way clusters. We did not 
explore a dependence of binding energy on projected galacto- 
centric radius due to distance uncertainties. 

6. CONNECTION TO OTHER SPHEROIDAL SYSTEMS 

Many efforts over the years have investigated the connec- 
tion between globular clusters and other spheroidal systems. 
O ne valuable approach is the /c-space formalism developed 
hv lBenderet"atl 17992). This space is comprised of three or- 
thogonal axes that are proportional to mass (k\), mass-to-light 
ratio (kt), and a third perpendicular axis that sc ales as surface 
brightn ess cubed times mass-to-light ratio (K3). B urstein et alJ 
( 1997) compiled K-space coordinates for a large number of 
spheroidal systems, with masses from globular cluster scales 
to clusters of galaxies. For globular clusters, the K-space co- 
ordinates are derived as 

Ki = (log a 2 + log r,)/ V2 + 0. 1 1 (4) 
k 2 = (log a 2 e + 2 log/, -log r e )/V6 + 0.06 (5) 

and _ 

k 3 = (log a 2 -log/,- log r e )/V3 + 0.09, (6) 

where a e s» cr^, r e « r/,, and I e is the mean surface brightness 
within r e in units of fi-band solar luminosities per square par- 
sec. 

In Figure [9] we plot the three K-space projections of the 
fundamental plane for NGC 5128 clusters (open circles), 
along with Galactic globular clusters (crosses), those in M3 1 
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Fig. 7. — Core and half-light fundamental plane correlations for globular clusters. The core fundamental plane projections (top panels) 
demonstrate that the NGC 5128 clusters (open circles), including those with possible extratidal light (partly filled circles), approximately 
follow the same relations as globular clusters in the Local Group. These globular cluster systems include those for our Galaxy (crosses), M31 
(stars), M33 (squares), the Magellanic Clouds (triangles), and the Fornax dwarf galaxy (pentagons). Half-light parameters (bottom panels) 
show comparable relations. Data were culled from the literature sources cited in the text. 



(stars), and galaxies (diamonds and open triangles, data from 
lBursteinetalll997l) . This figure indicates that the NGC 5128 
clusters are quite similar, although systematically more mas- 
sive, than the Galactic clusters. In the K\ coordinate, globular 
clusters are clearly separated from galaxies. This does not 
necessarily rule out the existence of objects in this gap, as a 
simple extrapolation of the mass functions for both globular 
cluster systems and dwarf galaxies suggest that such objects 
should be extremely rare and, in the case of dwarf galaxies, 
quite difficult to detect. Simple scaling arguments suggest 
that globular cluster systems in massive ellipticals such as 
NGC 5128 are arguably the best place to identify objects in 
this region, provided there is not some other physical mecha- 
nism that sets a maximum globular cluster mass. The globular 
clusters in the present study are a step toward eliminating the 
mass gap. 



There is also a separation in the K3 coordinate (M/L) be- 
tween globular clusters and galaxies, although not as stark 
as n\ (see Figure |9] left panel). The mean mass-t o-light 
ratio of Galactic gl obulars is 1.45 (lMcLaughlinll2000h . and 
Larse n etaLl (|2002) found a comparable number for M33. 
While the higher mass-to-light ratio for NGC 5128 clusters 
is somewhat surprising, this may be a common feature of 
more massive clusters. For example, the mass-to-light ratios 
of wCen and Gl are in the same range as those in NGC 5128, 
as are the mass-to-light ratios (or K3 values) of other massive 
clusters in the Galaxy and M3 1 . The K3 values for dwarf and 
regular elliptical galaxies are higher than for most globular 
clusters, reflecting their larger mass-to-light ratios, although 
they overlap with the most massive globular clusters. The fun- 
damental plane relations hip for spheroidal galaxies re quires 
T v oc Z°- 2 Jvan der Marell199lHMagorrian et alii 9 98). A re- 
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7. DISCUSSION 
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Fig. 8. — Binding energy vs. absolute V magnitude for globular 
clusters. The NGC 5128 globular clusters (circles) fol low approx- 
imately the same relation found by McLaughlin (2000) for Galac- 
tic clusters (crosses), although they diverge somewhat toward higher 
binding energy for a given My. Clusters from M31 (stars), M33 
(squares), and the Fornax dwarf galaxy (pentagons) are also shown. 
u;Cen and M31 Gl are shown as an emphasized cross and star, 
respectively. NGC 5128 clusters with possible extratidal light are 
marked with partly filled circles. 



lationship between mass (or luminosity) and mass-to-light ra- 
tio for the most massive globular clusters is suggested by their 
apparently higher K3 values. 

The K2 (Ig x M/L) coordinate is essentially an indicator of 
surface brightness at the effective radius as its dependence on 
M/L is much weaker than I e , K2 exhibits a great deal of over- 
lap in the properties of all spheroidal systems. In fact, the 
galaxies overlap much more closely with Galactic globular 
clusters than those in NGC 5128. This is because K2 is dom- 
inated by the mean surface brightness within the effective ra- 
dius and the NGC 5128 clusters are on average lower in mean 
surface brightness within this radius (e.g. see Figure[7}. 

While the NGC 5128 globular clusters overlap more with 
other globular clusters than galaxies, they are an even bet- 
ter match t o the nuclei of the dwarf ellipticals studied by 
iGeha et all (120021) . who measured both integrated and nu- 
clear velocity dispersions and surface brightnesses for a sam- 
ple of Virgo dwarf ellipticals. While the integrated properties 
of nucleated dwarf galaxies (large, open squares) fall within 
the region of K-space inhabited by other dwarf ellipticals 
and spheroids (diamonds), their nuclei (squares with crosses) 
overlap very well with the NGC 5128 clusters. 

The H1 — H2 relation for giant elliptical galaxies and bulges 
exhibits an anticorrelation between mass and surface bright- 
ness that sets them nearly orthogonal to the dwarf galaxies, 
which are correlated in K\ - «2- This correlation reflects the 
correlation between lumi nosity and larger co res of higher sur- 
face brightness found by Kormendv ( 1985). While this cor- 
relation suggests mass loss due to winds, Bender et a 
found that mass loss would produce too steep a relation be- 
tween K\ and K2, instead suggesting the observed dwarfs are 
the remnants of a larger, unobserved population. The massive 
globular clusters appear to exhibit a similar K1 — K2 correla- 



The strong dynamical and photometric similarities of these 
massive globular clusters to the nuclei of dwarf galaxies has 
interesting implications for models which posit stripped dwarf 
galaxy nuclei as the origin of the most massive globular clus- 
ters. The most massive Local Group clusters, such as wCen 
in our own Galaxy and Gl in M31, have both been dis- 
cussed and modeled as stripped dwarf nuc lei (Freeman 1993; 
[ Mevlan et afl 1200 ll iGnedin et all 120021 iBek ki & Free man! 
2003). The position of wCen and Gl in Figures 0and|9]show 
that the NGC 5128 clusters occupy a similar region of the 
fundamental plane and K-space. All of these massive clus- 
ters have comparable masses, mass-to-light ratios, and central 
surface brightnesses. 

Support of the interpretation of these clusters as tidally- 
stripped galaxy nucl ei is provided by th e tentative detection 
of extratidal light by Harris et al. (2002) for six of the globu- 
lar clusters (marked with an x in Table |2] and [5] partly filled 
circles in the figures). In the context of the stripped-dwarf 
model, it is tempting to view this extratidal light as the last 
vestiges of the extended dwarf envelope around these nuclei, 
although extratidal light may also be present due to the evap- 
oration of cluster stars by two-body relaxation. Deeper ob- 
servations of these and other massive clusters would be ex- 
tremely valuable to confirm and quantify this extratidal light. 
Such observations would also better quantify the ellipticities 
of these globular clusters, which are comparable to those of 
wCen and Gl but greater than those of dwarf nuclei. 

These globular clusters may also simply represent the up- 
per end of the globular cluster mass function in NGC 5128. 
The similarities between these globular clusters and the most 
massive globular clusters in the Local Group are more es- 
tablished than the interpretation of the better-studied Local 
Group clusters as tidally-stripped dwarf nuclei. In addition, 
there is evidence for yet more massive young star clusters that 
demonstrate that there is overlap between the masses of the 
most massive star clusters and the least massive dwarf galax- 
ies. The most massive star cluster known is W3, a clus- 
ter with a = 45 km s " 1 in the merger remnant NGC 7252 
( iMaraston et al J 120041) . The inferred mass of this cluster is 
nearly 10 8 M©, even m ore massive than the nuclei of the 
dwarf galaxies studied by Geh a et alJ ( 120 02). Even if we were 
to make the extreme suggestion that all of the NGC 5 128 glob- 
ular clusters are in fact relic dwarf nuclei, the existence of W3 
demonstrates there is overlap between the masses of galaxies 
and star clusters. 

On the other side of the mass gap, the ultra- compact dwarf 
galaxies in the Fornax cluster (Drinkwa ter et al.l20 03) are ac- 
tually less massive than W3 and overlap with the most mas- 
sive NGC 5128 clusters. These dwarfs have velocity disper- 
sions of cr = 24-37 km s" 1 , effective radii of 10 - 30 pc, 
masses in the ra nge 10 7 ~ 8 and mass -to-light ratios of 2 - 
4 in solar units ( iDrinkwater et alJl2003l) . They are thus very 
comparable to the larger and more massive of the NGC 5128 
globular clusters. Numerical models for these ultra-compact 
dwarfs have shown that they can form from nucleated dwarfs 
that ha ve been stripped of their envelopes by tidal forces in a 
cluster ( Bekki et al. 2003). This stripping process can explain 
the ultra-compact dwarfs in Fornax, although it will operate 
over approximately a factor of 2 smaller radius in a smaller 
group like NGC 5128 due to decreased tidal shear. The radius 
for the tidal stripping of NGC 5128 globular cluster progeni- 
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Fig. 9. — Projections of 3-dimensional K-space. The NGC 5128 globular clusters (open, large circles, or partly filled circles for those with 
possible extratidal light) are in the region of K-space that overlaps the most massive Galactic globular clusters (crosses), M3 1 clusters (stars) 
and dE nuclei (squares with crosses). ujCen and M31 Gl are shown as an emphasized cross and star, respectively. The integrated properties of 
the dEs (open squares) place them among other dEs and dSphs (open triangles). Elliptical galaxies and bulges are also shown (diamonds). 



tors is likely to be on order 10 kpc ( Bek ki et al.l2003l see their 
Figure 7), comparable to the projected distances of 5 - 23 kpc 
for some of the NGC 5128 globular clusters in this study. It is 
therefore plausible that tidal stripping could have transformed 
at least some nucleated dwarf galaxies into these NGC 5128 
globular clusters. 

Another aspect of the tidally-stripped dwarf model is that 
their dark matter halos cannot be too cuspy as otherwise the 
relatively concentrated dark matter core will be too effective 



at reta ining the stellar envelope. In their study, BekkLetaD 
d2003|) used the dark matter p rofile of ISalucc i & Burkert 
(2000) (originally proposed by Burkert (1995) for dwarf 
ga laxies), rather than the more co mmonly adopted profile 
of iNavarro. Frenk. & White! (119961) . because it has a flatter 
core. However, this requirement for a relatively flat dark mat- 
ter core may remove one possible explanation for the higher 
mass-to-light ratio of these massive globular clusters. While 
a residual dark matter halo from their dwarf galaxy past is 
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a possible explanation for the larger mass-to-light ratios of 
these NGC 5128 clusters, a dark matter profile with a flat, 
low-density core w ill also be more efficiently stripped away 
( BekJ d~etalll2003|) . An alternative explanation for the high 
mass-to-light ratios of these massive globular clusters is if 
they for med in a starburst wi th a truncated stellar initial mass 
function ( Chariot et a 1. 1993). A mass function with relatively 
more low-mass stars would produce a higher mass-to-light ra- 
tio. 

The detailed shapes and kinematics of the nuclei of nucle- 
ated dwarfs and the most massive globular clusters may pro- 
vide one way to further investigate the potential connection 
between these two populations. T he most massive g lobular 
clusters have significant ellipticities (Harris et al. 2002), while 
this do es not appear to be the case for the nuclei of nucleated 
dwarfs (Geha et al. 2002). However, the same tidal forces that 
strip a dwarf envelope may also induce significant elliptici- 
ties. The importance of rotational flattening and anisotropies 
in the velocity distribution may also serve to distinguish be- 
tween dwarf nuclei and globular clusters. 

The hypothesis that some of the most massive globular clus- 
ters are the nuclei of galaxies offers an appealing explanation 
for recent evidence of an inte rmediate-mass black hole in Gl 
(Gebh ardt. Rich. & Hoi 12002). The mass of this black hole 
falls on the same Mbh - c relationship for galaxies and sug- 
gests that the formation mechanisms for black holes in star 
clusters and galaxy spheroids are similar. If Gl is instead sim- 
ply a tidally-stripped, nucleated dwarf galaxy, the problem is 
reduced in complexity and only one physical mechanism for 
black hole growth is required to explain the Mbh - & relation. 

An alternate interpretation of the similarity between the 
most massive globular clusters and the nuclei of nucleated 
dwarfs is that the later are simply star clusters that have mi- 
grated to or formed at the centers of these dwarfs. The prop- 
erties of these objects would then be probes of massive star 
clusters in different environments, rather than of actual over- 
lap between the properties of the most massive star clusters 
and the least massive galaxies. This interpretation also ex- 
plains their similar location in K-space, although stands in 
contrast to the simple explanation for the intermediate-mass 
black hole in Gl. 

8. SUMMARY 

We have measured velocity dispersions for a sample of 
14 globular clusters in the nearby, luminous elliptical galaxy 
NGC 5128, the first such study of the globular cluster system 
of a luminous elliptical galaxy and the first such study outside 
the Local Group. These clusters have velocity dispersions in 
the range 10 - 30 km s" 1 , comparable to the largest previ- 
ously measured values for globular clusters. We have used 



measured King model structural parameters for these clusters 
from the literature to derive masses for all 14 clusters. These 
clusters are comparable in mass to the most massive Galactic 
globular cluster uiCen and M3 l's Gl . From these data we find 
the following: 

1. The globular clusters in NGC 5128 approximately fol- 
low the same fundamental plane relationships as Local Group 
globular clusters and extend them to approximately an order 
of magnitude higher mass and luminosity. 

2. The mean mass-to-light ratio of these clusters is larger than 
for typical Local Group globular clusters, although compara- 
ble to the more massive Local Group clusters. 

3. These clusters begin to bridge the mass gap between the 
most massive globular clusters and the least massive dwarf 
galaxies. In particular, there is very good overlap in the pho- 
tometric, structural, and kinematic properties of these clusters 
and the properties of both nucleated dwarf elliptical nuclei 
and ultra-compact dwarf galaxies. 

4. The large masses of these clusters, combined with the 
possible detection of extratidal light for some objects by 
lHarris et all J2002I) . suggest that some of these globular clus- 
ters are in fact the nuclei of tidally stripped dwarf galaxies. 

The common properties of the most massive star clusters 
and the nuclei of the least massive dwarfs suggest that both 
the formation mechanisms for star cluster and galaxies can 
produce objects in the same region of the fundamental plane 
or K-space. Alternately, the nuclei of nucleated dwarf galax- 
ies may simply be star clusters that happen to lie in the cen- 
ters of galaxies, rather than true galaxy nuclei. Future spec- 
troscopic observations of additional massive globular clusters 
could quantify the relative contribution of relic dwarf nuclei 
to this population through kinematics and with stellar popula- 
tion models, while deep, high-resolution images could pro- 
vide better measurements of structural parameters, particu- 
larly in the core, and search for and study diffuse, low surface- 
brightness envelopes. 
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